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Abstract. We prove finite-time blowup for spherically symmetric and nega- 
' tive energy solutions of Hartree-Fock and Hartree-Fock-Bogoliubov type equa- 

I tions, which describe the evolution of attractive fermionic systems (e. g. white 

f™*^ , dwarfs). Our main results are twofold: First, we extend the recent blowup 

^ . result of [Hainzl and Schlein, Comm. Math. Phys. 287 (2009), 705-714] to 

' Hartree-Fock equations with infinite rank solutions and a general class of New- 

tonian type interactions. Second, we show the existence of finite-time blowup 
for spherically symmetric solutions of a Hartree-Fock-Bogoliubov model, where 
an angular momentum cutoff is introduced. We also explain the key difficulties 
' encountered in the full Hartree-Fock-Bogoliubov theory. 
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1. Introduction 



< 

1-^^ ' In this paper, we study time-dependent generalized Hartree-Fock equations that 

, arise as quantum fermionic models for the dynamics of relativistic stars (like white 

dwarfs and neutron stars). 

Here, our particular interest is devoted to the existence of finite-time blowup 
solutions, which physically indicate the onset of "gravitational collapse" of a very 
massive star due to its own gravity. Indeed, such a catastrophic scenario was already 
f^f^ • predicted by the physicist S. Chandrasekhar in 1930's, based on an intriguing and 

, heuristic combination of special relativity and Newtonian gravity. However, it is fair 

' to say that a rigorous understanding of the dynamics of the gravitational collapse 

\ of a star - in terms of relativistic quantum mechanics and Newtonian gravity to 

0^ ' start with - is still at its beginning. In fact, the proposed mathematical models 

\ (such as the generalized Hartree-Fock equations considered below) are formulated 

' as nonlinear dispersive evolution equations, which exhibit the essential feature L^- 

mass-criticality and thus reflecting the physical fact of a so-called Chandrasekhar 
^ limit mass 

■ A/* ~ 1.4Msu„, 

■ beyond which gravitational collapse of a cold star can occur. For a detailed discus- 
sion of Hartree-Fock type theory in astrophysics, we refer to [inill^] and references 
given there. See also [111 HT] for rigorous work on Chandrasekhar 's theory in the 
time-independent setting. 

Let us now briefly outline the main content of the present paper. As a relativistic 
quantum model of the dynamics of star, we consider time-dependent generalized 
Hartree-Fock equations with relativistic kinetic energy and an attractive two-body 
interaction of Newtonian type. More specifically, we study equations of Hartree- 
Fock (HF) and Hartree-Fock-Bogoliubov (HFB) type, and our main results on finite- 
time blowup of spherically symmetric solutions for these models will be stated in 
Sections [5] and O respectively. 
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In fact, this paper aims at two directions: First, we consider the HF equations 
and we extend the recent blowup results of 12 to solutions with infinite-rank den- 
sity matrices and a broader class of two-body interactions. Second, we address the 
finite-time blowup for HFB type evolution equations, which also incorporate the 
physically important phenomenon of Cooper pairing for attractive fermion systems. 
However, it turns out that the analysis of HFB equations brings in new key difficul- 
ties, which were completely absent for HF. Most importantly, the nonconservation 
of higher moments of angular momenta in HFB theory makes an adaptation of the 
strategies used in |10[ I12j an insurmountable task at the present. To circumvent 
this difficulty, we introduce an HFB model with an angular momentum cutoff and 
prove a finite-time blowup result for this simplified model; see Section [3] for more 
details. In the long run, we do hope that our arguments developed in this paper 
will be of use to understand the finite-time blowup for the full-fledged HFB model 
without any cutoffs. 

Outline of the paper. This paper is organized as follows. In Sections [2] and 
131 we first collect some preliminary facts and we then state our main results on 
finite-time blowup for spherically symmetric solutions of HF and HFB type models. 
Furthermore, we explain the key difficulties encountered for HFB theory. The proofs 
of the main theorem will be given in Section [H and Appendix A contains the proof 
of a technical lemma. 

Conventions and notations. For the physically inclined reader, we remark that 
we employ units such that Planck's constant h and the speed of light c both satisfy 
h = c = 1. As usual, the letter C denotes a positive constant depending only on 
fixed quantities (such as initial data etc.) and C is allowed to change from line to 
line. 
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2. Blowup in Hartree-Fock theory 

In this section, we consider the Hartree-Fock (HF) model describing a system 
of pseudo-relativistic fermions with average particle number A > and attractive 
interactions of Newtonian- type (like classical gravity). With regard to physical 
applications, we may think of such system representing a white dwarf star, where 
effects due to special relativity play a significant role but gravity can still be treated 
in the Newtonian framework. Of course such a model lacks Poincare covariance, 
and hence it is usually referred to as pseudo-relativistic. Nonetheless, many essen- 
tial effects (such as the existence of a Chandrasekhar limit mass of white dwarfs) 
are already predicted - even quantitatively within good approximation - by such 
pseudo-relativistic models. Indeed, it should not go unmentioned here that S. Chan- 
drasekhar 's original and acclaimed physical theory of gravitational collapse rests on 
a pseudo-relativistic model of a star in its semi-classical regime; see [?]• 
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According to HF theory, the state of the system is described in terms of a density 
matrix 7 which is a self-adjoint operator acting on satisfying ^ 7 ^ 1 

(reflecting the Pauh principle for fermions) and the normalization condition Tr(7) = 
A fixing the (average) number of particles. 

Let V denote the interaction potential between the particles and, for notational 
convenience, let k > be a coupling constant controlling the strength of the inter- 
action. Then the corresponding HF energy reads 

(2.1) £HF(7) = Tr(X7) + | JJ V{x~y)(^p^{x)p,{y)-\j{x,y)\'^dxdy. 

Here 7(2;, y) is the integral kernel of the trace-class operator 7 and p-yix) = 7(2;, x) 
is the associated density function. The pseudo-differential operator 

K = V-A -1- to2 

describes the kinetic energy of a relativistic quantum particle with rest mass m ^ 0. 

In the present section, we are interested in the corresponding time-dependent 
HF equations, which can be formulated as the following initial-value problem: 



(2.2) 

Here [A, B] = AB — BA denotes the commutator and 




(2.3) := + m? + k{V * pj) - kV{x - y)^{x, y) 

is the so-called mean-field operator which depends on 7 and acts on the one-body 
space L^(M'^;C). Here and henceforth, the symbol * stands for convolution of 
functions on 'M? . For expositional convenience, we use a slight abuse of notation 
by writing V{x — y)j{x,y) for the operator whose integral kernel is the function 
{x,y) I— *■ V{x — y)^{x,y). We remark that the number of particles Tr(7t) and the 
energy EnFi'^t) are both conserved along the flow given by (j2.2p . as we will detail 
later on. Also, the appropriate set /Chf of initial data for the evolution equation 
(H^ will be defined below. 

More specifically, we consider radial interaction potentials V = V{\x\) that are 
essentially Newtonian at small and large distances. By this, we mean that V is of 
the form 



(2.4) 









^w{\x\). 







Here w is a smooth and decaying function, satisfying the following conditions that 
we impose throughout this paper: 

(2.5) w£C^[0,oo), r|w(r)KC and {1 + r^)^+'' \w' {r)\ C 

for all r ^ and some e > 0. Physically, one may think of w as describing screening 
effects due to internal degrees of freedom of the gravitating fermions. Indeed, as 
remarked above, our arguments carry over to the case where spin is taken into 
account, and w = w(|a;|) would he a q x q matrix instead of a multiple of the 
identity (with an adequate redefinition of the HF-energy function £hf (7) above) . 

In the special case of 7 being an orthogonal projection of finite rank (i. e, we have 
7 = X]i=i {''Pi) {''PiD and vanishing w = 0, the existence of finite-time blowup for the 



time-dependent HF equation (|2.2p was recently studied by Hainzl and Schlein in 



"'^For simplicity, we discard the spin of the particles throughout. But all our arguments can 
be easily generalized to particles having q internal degrees of freedom, where L^(K''; C) has to be 
replaced by L'^(R^;C'') etc. 
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[T^ , by extending the arguments developed by Frohhch and Lenzmann in [TU] . The 
main result of 12J established that any (sufficiently regular and rapidly decaying) 
initial datum 7 with spherical symmetry (see Definition [1]) and negative energy 
^hf(7) < leads to finite-time blowup of the solution 74. Note that the assumption 
of finite rank was important for the method used in |12] to work. In this paper, 
we present a new technical approach which allows us to prove blowup for 7 having 
infinite rank. Moreover, we show here that the existence of blowup solution is 
stable with respect to a certain class of perturbations w of the Newton interaction, 
satisfying (|2.5p . 

After this preliminary remarks, we now state our results on the HF equation 
p.2p . To this end, we denote by &p the Schatten space (see [201 119) of operators 
A acting on L'^{R^,C) such that ||^||^ = Tt\A\p < 00. Further we introduce a set 
of fermionic density matrices 

(2.6) /Chf := {7 = 7* G -^hf : ^ 7 1} , 
where the Sobolev-type space Ahf is defined by 

(2.7) A-HF {7 e 61 : hlUnp < f»} 
with the norm 

(2.8) I|7|Uhp:= (1-A)1/S(1-A)1/4 . 

61 

In fact, it turns out that the initial- value problem (j2.2p is locally well-posed in /Chf. 

Theorem 1 (Well-posedncss in Hartrec-Fock theory). Suppose that w satisfies 
p.Sp . For each initial datum 70 G /ChFj there exists a unique solution 74 G 
C°([0, T), /Chf) H C^([0, T); A'^p) with maximal time of existence < T ^ co. 
Moreover, we have conservation of energy and expected particle number: 

'^^HF(7t) = S-apilo) and Tr(7t) = Tr(7o) for sC < < T. 

Finally, if the maximal time satisfies T < 00, then Tr(— A)"'^/^7( -^00 as t ^ . 

This theorem was proved in [lOj in the special case w = and for 7 with finite 
rank. However, it is straightforward to extend the arguments of [TD] to the case 
w 7^ and 7 with infinite rank; see also [HIS]. 

Remark 1. i) By using the conservation laws, we can deduce that sufficiently small 
solutions 7t extend to all time. More precisely, there exists a universal constant 
C* > such that 

Tr(7o)K3/2 < 

implies global- in-time existence and a priori bounds in energy space; i.e., we have 
T = +00 and that ||7t||A'HF ^ for all times t ^ 0; see [T3]. Physically, this global 
wellposedness result states that white dwarfs of sufficiently small mass have a well- 
defined global dynamics and, in particular, no gravitational collapse can occur. 

ii) Particular global-in-time solutions of (|2.2p are stationary states satisfying 
[Hj,^] — 0. Important examples for stationary states are given by the minimizers 
of She (7), subject to the constraint Tr(7) — A with Xtc'l^ not too large, which were 
proven to exist in [Tl] . 

iii) Well-posedness can also be shown to hold true in Sobolev-type spaces 'W 
with the norm ||7||w= = ||(1 - - A)''/2||g^ for s ^ 1/2, provided the ini- 
tial datum 7o belongs to 'W . Indeed, we will make tacitly use of this fact about 
this persistence of higher regularity, guaranteeing the well-definedness of our calcu- 
lations below. 

We now define spherically symmetric states for which we will be able to prove 
blowup under some assumptions: 
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Definition 1 (Spherically symmetric HF states). We say that 7 G /Chf is spheri- 
cally symmetric when ^(Rx, Ry) = 7(x, y) for all x,y and all R G 50(3). 

It is easy to verify (see Lemma U below) that spherically symmetry of 74 is 
preserved under the flow (|2.2p . We also note that, if 70 is sufficiently regular, 
then the condition of spherical symmetry can also be written as the commutator 
condition 

(2.9) [luL]^Q. 

Here L = —ix A V^: is the angular momentum operator, and A denotes the cross 
product on R"^. 

Our main result of this section is the following blowup criterion. The proof will 
be given in Section below. 

Theorem 2 (Blowup in Hartree-Fock theory) . Suppose that w is a radial Junction 
satisfying (|2.5p and that k > 0. Let 70 G /Chf &e spherically symmetric and suppose 
that 

Tr |x|So + Tr (-A)7o + Tr \L\^-fo < 00 

where L — ~ix A V^; denotes the angular momentum operator. 

Our conclusion is the following: If 70 has sufficiently negative energy, that is 

(2.10) SMlo) < - J (Tr(7o))' sup\w{r) + r w' {r)\_ 

(where |/|_ ^0 denotes the negative part of f), then the corresponding solution 74 
to S2. 2\) blows up in finite time; i. e., we have Tr(— A)^/^7t ^00 as t ^ T~ 
for some T < 00. 

Remark 2. By scaling, it is simple to show the existence of 70 G /Chf satisfying 
the condition ()2.10|) for sufficiently large coupling constants k, or equivalently for 
a sufficiently high particle number Tr(7o). Note also that Tr(7o)K'^/^ must be 
sufficiently large, since otherwise T = +00 holds, as mentioned in Remark [T] above. 

Regarding the proof of Theorem [21 we remark that we follow a virial-type ar- 
gument developed in [TO] for the Hartree equation, which was recently extended 
in [T2] to Hartree-Fock equations, by using an additional conversation law for the 
square of the angular momentum Tr(|Lp7t) for radial solutions. More specifically, 
it was shown in 12J that spherically symmetric solutions 74 of the Hartree-Fock 
equation (12. 2p with V = —l/\x\ satisfy the inequality 

(2.11) Tr(M7t) < 2£hf(7o) + CTr(7o)Tr((l + \L\''-fo)t + C, 

provided that 70 is finite-rank. Since the relativistic "virial" operator 

3 

M = ^ Xi\/-A + rn^Xi 

is nonnegative, this yields a finite maximal time of existence T < 00 and hence the 
desired blowup result. In the arguments developed below, we remove the finite-rank 
condition, by expanding the solution 7^ in terms of sectors of angular momentum. 
This approach is also essential when addressing the blowup for HFB theory, as done 
in the next section. However, the conservation law for Tr(|ip7t), which plays an 
essential role in the argument sketched above, is no longer at our disposal in HFB 
theory, as will be detailed below. 
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3. Blowup in Hartree-Fock-Bogoliubov theory 

In this section, we now turn to the Hartree-Fock-Bogohubov (HFB) model, which 
generahzes the well-known Hartree-Fock model studied in the previous section. It is 
fair to say that HFB theory is a widely used tool for (partially) attractive quantum 
systems, which are found - for example - in nuclear and stellar physics; see [TOIIHIIB] 
for a detailed exposition. As a notable fact, we mention that HFB yields (in some 
limiting approximation) the well-known Bardecn-Cooper-Schrieffer (BCS) theory, 
which plays a fundamental role for the understanding of superconductivity and 
superfluidity. 

The outline of this section is briefly as follows. First, we introduce the HFB 
energy and its associated time-dependent equations. Second, we point out the key 
difficulties in the blowup analysis encountered within HFB theory. Finally, we 
state our main result on finite-time blowup for the HFB model with an angular 
momentum cutoff; see Theorem |4] below. 

3.1. HFB Theory: Preliminaries and Difficulties. 

3.1.1. HFB Energy. The HFB model involves two variables: the same one-body 
density matrix 7 as in HF theory, and an additional two-body fermionic wavefunc- 
tion a{x,y) £ L^(M'^;C) A L^(Ili'^;Cjl describing the so-called Cooper pairs. The 
amount of Cooper pairing is given by // \a{x,y)\'^dx dy; and it must vanish when 
7^ = 7 is a pure Hartree-Fock state. For simplicity, we will always use the same 
notation for the two-body wavefunction a and the corresponding (non self-adjoint 
but compact) operator acting on the one-body space L^(M^; C), whose integral ker- 
nel is the function {x,y) 1— > a{x,y). Note that the assumption that a is fermionic 
can then be expressed as = —a. We remind the reader that again we have 
neglected any spin degrees of freedom to simplify our presentation. 

Within the HFB framework, the operators 7 and a are, by assumption, linked 
through the following operator constraint: 

(3.1) ""^^{a* 1^7)^ (J ?) onL^(M3;C)eL^(R3;C). 

In fact, this inequality guarantees that the pair (7, a) is associated with a unique 
quasi-free state in Fock space; see [2]. The average number of particles is given (as 
in HF theory) by the trace Tr(7). 

The HFB energy functional, for a system interacting through a radial two-body 
potential V, is given by 

(3.2) £hfb(7,«) :=^hf(7) + 5 // \aix,y)\''Vi\x - y\) dx dy. 

It is then obvious that a ground state (if it exists) always has a = in the case 
of a purely repulsive potential V ^ 0. Therefore, the use of HFB theory is only 
meaningful when V is not everywhere positive. In fact, in this case, it is often 
observed that there is a phase transition from an HF ground state to an HFB 
ground state when the coupling constant k is increased. 

As in the previous section, we consider radial potentials V — V{\x\) of the form 
()2.4p above. Furthermore, we remark that the HFB energy £HFB(7,a) can be 
appropriately defined (see [141 ) on the set of generalized density matrices: 

(3.3) /Cheb = { (7, a) £ ^heb : ([j [j) ^ ( J. 1 ! ^) ^ ( J j) } ■ 



'Here A denotes the antisymmetric tensor product. 
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The space Ahfb is the correct energy space for generahzed density matrices, which 
is defined as 

(3.4) Xhfb = {{l,a) <E &i X &2 : 7* = 7, = -a, \\{-f,a)\\x-tiFB < 
and is equipped with the norm 

(3.5) ||(7,a)|UHPB = 11(1 - A)i/S(l - A)'/le. + 11(1- ^)'^'"lle.- 

We remind the reader that ©1 and 62 denote the space of trace-class and Hilbert- 
Schmidt operators on L^(R^;C), respectively. Also, note that ||(1 — A)^/'*a||e2 is 
equivalent to the Sobolev norm ||a(-, •)||^i/2(^3xr3) where a{x,y) is regarded as a 
two-body wavefunction defined on R'^ x M^. 

The existence of minimizers for the HFB energy fHFB(7, ct), under the constraint 
that Tr7 = A for Ak^^^ not too large, was given in [H]. Moreover, some important 
properties were derived in [TJ I14j . We note it is not known whether a ground state 
always has a 7^ for V taking the form ()2.4p , but this is expected to hold true (at 
least for sufficiently large local potentials w). 

In view of p.ip and following [HIH], we introduce a generahzed one-body density 
matrix 

(^•«) 1-7)' 
as well as the corresponding HFB mean-field operator 

which acts on L^(R^; C) ® L^{W^; C), with denoting the conjugate of H^. Here 

(3.8) -.^ K + k{V * p.y){x) - KR^{x,y) 

denotes the Hartree-Fock mean-field operator introduced before in ()2.3p ; and Hq, is 
given by its kernel 

^a(x, y) = kV{x - y)a{x, y). 
Moreover, it turns out to be convenient to define the number operator as 

(3.9) AA^(J 

Note that F commutes with M if and only if a = 0, i. e., if and only if the corre- 
sponding quasi-free state in Fock space also commutes with the number operator 
A/"; see, e. g., [2]. 

Remark 3. It was shown in [21 |T3] that any minimizer of the HFB energy p.2p 
solves the following self-consistent equation: 

(3.10) T = X(-oofl){FT-^l^^)+D, 

for some Lagrange multiplier /i < chosen to ensure the constraint Tr(7) — A, 
and where _D is a finite rank operator of the same matrix form as F, satisfying 
ran(£') C ker(Fr — yJ^)- It was also shown in O [T] that if a 7^ and w = 0, then 
7 must be infinite rank. 
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3.1.2. The HFB evolution equation. Let us now turn to the time-dependent HFB 
theory. This evolution equation can be concisely written in terms of the generalized 
density matrix F — r(7, a) introduced in ()3.6|1 . The corresponding initial-value 
problem reads as follows: 

(3.11) \h¥^=^^^ 

[ (7|t=o,a|t=o) = (7o,ao) e ^HFB- 

This evolution equation can also be regarded as obtained from "forcing" the many- 
body Schrodinger evolution in Fock space to stay on the manifold of HFB states, 
by means of the Dirac-Frenkel principle in quantum mechanics. With regard to the 
wellposedness of its initial- value problem, we note that a straightforward adaptation 
of [HI [51 m [ini [m [in yields the existence and uniqueness of a maximal solution. 
We record the wellposedness fact as follows. 

Theorem 3 (Well-posedness in Hartree-Fock-Bogoliubov Theory). Assume that w 
is a radial function satisfying (|2.5p and that ^ k < A/tt holds. For each initial 
datum (70, ao) G A^hfB; there exists a unique solution (74, at) G C''([0, T); /Chfb) H 
C"'^([0, T); A'^pg) solving (|3.1ip with maximal time of existence < T ^ oo. More- 
over, we have conservation of energy and expected number of particles, i. e., 

Silt, at) = £(70, ao) and Tr(7f ) = Tr(7o) for ^ t < T, 

as well as the a priori bound for the pairing wave function: 

TT{a;at) C for 0^t<T. 

Finally, if T < 00 holds, then we have Tr(— A)"'^/^7t —> 00 as t ^ T~ . 

Remarks 4. i) It was shown in [13] that when Tr(7o)K'^/^ is sufficiently small, then 
one has T — +00 and (74, at) is bounded in A'hfb- 

ii) Note the condition k < 4/7r as opposed to HF theory; see Section [3.1.31 for 
more details on this. 

iii) It is straightforward to verify that if the generalized matrix Fq associated with 
the initial condition (70, ao) is an orthogonal projector on L^(R'^;C) ® L^(R"^;C), 
then so is Fj for all times t e [0,T). Similarly, if ao (i.e., the initial datum 
is a Hartree-Fock state) then we have at = for all time. Therefore the equation 
p. lip is a generalization of the HF evolution equation (|2.2p . 

iv) The a priori control on Tr (a j at) simply follows from the operator inequality 
7^ -|- a*a ^ 7 (since F^ ^ F) and the conservation of Tr(7). See Section [3.1.31 for 
more details on the time evolution of a. 

v) As for stationary solutions of (j2.ip . we remark that it can be verified that if 
Fo satisfies [Fr^ -^A/", Fq] = for some n, then F(t) = e"'''^* Fq e*^^* is solution of 
equation p. lip . In particular, we see that HFB minimizers give rise to stationary 
solutions of p. lip , as one naturally expects. 

Remark 5 (Spin symmetry). There is a very natural spin symmetry in HFB theory, 
which was already used before in minimization problems for the Hubbard model in 
[5] and for general systems with an interaction potential of negative type in [T] . In 
the case of the Newtonian interaction, this was used in [H] . 

Let us assume only in this remark that the problem is now posed in L^(R^; C^) 
and let us introduce the unitary operators for k = 1,2,3 
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where ak are the Pauh matrices. A simple calculation shows that [Sfe,r] = is 
equivalent to [crfe, 7] = and akOL+ao}^, — 0. It can then be checkecHthat [E^, F] = 
for fc = 1, 2 if and only if = = 0, 7|t = 7|x, = a^i = and a|| = -a^i- 

Let us now assume that our initial state Fq takes the special form 
(3.12) 

^°={a^ "^''"^ 70 = TO® (J and ao = ao J 

where qq and tq act on L^(M''; C) with Tq = tq and = oq. If w is a multiple of the 
2x2 identity matrix, then both the kinetic energy K and our interaction potential 
do not depend on the spin. It is then clear that we have S^FrSfe ~ F^^T'S' for 
k = 1,2, 3. Hence we see that S^FtS^ S /Chfb solves the same time-dependent 
equation as Ft. By uniqueness we deduce that Ft has the same form as Fq for all 
times. Therefore if one seeks for solutions of the above special form, one can reduce 
the spin-1/2 case to a no-spin model with the condition = —a replaced by a 
singlet-type condition = a. All our results also hold true under this constraint. 
A similar reduction can be done for a system with q ^ 1 degrees of freedom. 

The special form (j3.12p is sometimes called time-reversal symmetry as the time 
reversal symmetry operators are usually defined through Tk = ^^'^k-^/^j^ _ Yi^K, 
where K is the conjugation (anti- linear) operator. Indeed, one verifies that TkTtT^^ 
solves the backward time-dependent equation. 

3.1.3. HFB: The Difficulties. In the previous section we have studied the blowup of 
solutions in the Hartree-Fock case. As we will explain now, the HFB generalization 
brings in new substantial difhculties as follows. 

• The pairing term coming from at is difficult to control. 

• Powers of the angular momentum Tr(|L|*7t) are generally not conserved 
along the HFB flow, even for spherical symmetric solutions (74, at) defined 
below. This is a striking contrast to HF theory. 

It is enlighting to rewrite the time-dependent equation (|3.1ip as a coupled sys- 
tem with a von Neumann-type equation for the one-body density matrix 7t and a 
two-body Schrodinger equation for the pairing wavefunction at- In fact, a simple 
calculation yields the following system: 



(3.13) 





r . d 

'Jt^' 


= [H^t : It] - 


- iGat I 




< 








+ KV{\x-y\) 








-'«((7t 


X + {■yt)y)V{\x ~ y\)^ at{x,y). 



The time-dependent equation for 74 contains the trace-class coupling term Ga aris- 
ing from the pairing wavefunction a. By an elementary calculation, this term is 
found to be Ga = iijlao:* — aH* ) or, in terms of its integral kernel, 



Ga{x, y) = i K / a{x, z)a{y, z) {Vljj — z) — V{x — z)) dz. 

Note that Tr(7f) is constant along the trajectories, because of Tr(GQ,) = 0. 

It is fair to say that the two-body wavefunction at evolves through a time- 
dependent two-body Schrodinger equation with Hamiltonian 

(3.14) (H,,)^ + {H,X + kV{\x - y\) - 4{jth + ijt)y)V{\x ~ y\) 



r for fc = 1,2,3 is indeed equivalent to 7|| = 7|| = 0, 7|| = 7|| 
and o = 0, as may be seen by a simple computation. We do not want to impose such a condition. 
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Here, the first two terms in this Hamiltonian form a one-body operator arising 
from the mean- field operator H^, which itself contains the potentials induced by 
7 as external sources. The third term KV^(|a; — y\) is a usual two-body Newtonian 
interaction term (perturbed by our local interaction w). The last term K,{^{''it)x + 
{"1t)y)V{\x — y\) in (I3.14|) is a non-self-adjoint (but compact) operator which is 
responsible for the fact that at does not have a constant norm along a given 
trajectory (except of course if at = 0). We note that since ^ k < 4/7r, the two- 
body operator ()3.14p is actually equal to a well-defined self-adjoint operator, plus 
a compact non-self-adjoint operator. Therefore, the two-body equation for at{x,y) 
is well-defined. 

Let us now detail the behavior of moments of angular momentum along the HFB 
flow. We recall for the reader's orientation that, in HF theory, higher moments of 
angular momentum are conserved, provided the initial state is spherically symmet- 
ric. This fact is one of the main tools used in blowup proof for HF, as we have 
already seen. First, we define spherically symmetric HFB states in a similar fashion 
as in the HF case. 

Definition 2 (Spherically symmetric HFB states). We say that (7, a) S /Chfb is 
spherically symmetric when ^{Rx, Ry) = 7(0;, y) and a{Rx, Ry) = a{x, y) for all 
eK^ q,^^ all R & SO{i). 

When the HFB initial datum (70,010) is spherically symmetric, it can easily 
be seen that (74,04) stays spherically symmetric on its time interval of existence. 
In terms of the angular momentum operator (provided that 7 and a are smooth 
enough) the condition of spherical symmetry may also be written as 

(3.15) [i,7t]=0 and [L., + Ly)at{x,y) = 0. 

Moreover, note also that we find that [|Lp,7t] ~ holds. However, we do not 
obtain any other information about the two-body wavefunction at- In particular, 
we do not have that {\La:\^ + \Ly\'^)at(x, y) — 0, which would mean that at{Rx, y) = 
at{x,Ry) = at{x,y) for all rotations R £ 5*0(3). 

Let us elaborate more explicitly the fact that Tr(|Lp7t) is not conserved for 
spherically symmetric solutions in HFB theory. To this end, we calculate the varia- 
tion of the expectation value of |Lp, assuming that [L, 74] = and {Lx+ Ly)at = 
for all times t G [0, T). We obtain 

(3.16) |Tr(|Lp7t) 



2 



l-^yP' 1 \-w{\x-y\) 



at 



As \x — y\^^ and w{\x — y\) do not commute with \Lx\'^ + l-^yPi the above term is 
usually not zero, and hence the square of the angular momentum is not conserved. 
This is a major difficulty that prevents us to use the same proof as in the HF case! 

To see that the square of the angular momentum does actually vary in time, 
let us consider a sufficiently smooth initial datum (70,00) with oq 7^ such that 
7o(a;,?/) = 7o(|a;|,|y|) and aQ{x,y) — aQ{\x\,\y\). The latter property means ex- 
actly that (I-LkP + \Ly\'^)aQ{x,y) = 0. It is clear from the above formula that the 
derivative of Tr(|Lp7t) vanishes at time i = 0. We have the following series of 
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equality: 

(|L,|2 + \Ly\^) ((70), + (70),) ao(|x|, \y\)V{x - y) 

= ((70). + (70),) {\L.? + \v\)V{x - v) 

= 2 ((7o):e + (7o),j) ia: • Lyao{\xl \y\)V{x - y) 
= 2 ((7oi)a; -Ly + L^- (7oi)j/) ao(kl, |2/I)^(a; - y) 
= 

since ^qL = Ljq = and (L^ + Ly)V{x — y) — 0. Using this, we can compute the 
second derivative of the expectation value of |Lp. We obtain 

(3.17) §,Tri\LUh., = .V^!!fN^, + \Ly^^^ ""d^l' I^D 

at^ ' \ \x — v\ 



^-y\ ' \X-y\ / L^^L^ 



+ K'{(3o{x,y),{\L,\^ + \Ly\^)(3o{x,y))^ 



where f3o{x,y) — ao(|a^|, |y|)u'(|a; — y\)- It is easy to see that p.l7p is always a 
positive quantity, except when = Q. 

To sum up: Even when we start at time t = Q with a state in the zero angu- 
lar momentum sector, it immediately acquires a change of the square of angular 
momentum, in stark contrast to HF theory 10 p. 743]. 

3.2. Blowup with Angular Momentum Cutoff". As we have explained, the 
second moment of the angular momentum, Tr(|Lp7t), is not conserved for spher- 
ical states in general. In fact, it seems quite difficult to get any useful bound on 
this quantity. There is yet another difficulty that we have not mentioned so far: 
the pairing term is difficult to control and our proof reveals that we actually need 
a bound on Tr(|L|^+^7t) for some £ > 0; see Lemma [7] below. It seems a very 
interesting question to investigate whether the angular momentum can indeed sub- 
stantially alter the collapsing behavior of a gravitational system in HFB theory, or 
actually prevent collapse at all (which, however, is very unlikely in our opinion). 

For all these reasons, we will now introduce a simpler model for which we are 
able to prove blowup of spherical solutions with sufhciently negative energy. The 
idea is to enforce an angular momentum cutoff in the model as follows. Let A ^ 
be a fixed integer, and define the (infinite-dimensional) orthogonal projectoi0 

^^A X[o,A(A+i)] (lip) , 

which projects ^^(M^; C) onto sectors with angular momentum £ ^ A(A + 1). We 
will now impose that our density matrices 7 and a satisfy for all times 

(3.18) Pa7 = 7^a=7, {PA(E)PA)a{x,y)^a{x,y). 

Note that the assumption on a is equivalent to {{Pa)x + {PA)y)oi{x, y) = 2a{x, y), 
or in operator terms to PaO£ = aPA = a. Our assumption p.l8p is the same as 
saying that we restrict ourselves to HFB states in Fock space !F that actually belong 
to the restricted Fock space 

^A = (Sn^q Af i3A C 

where S)a = PaL'^(^^', C). Hence, introducing the angular momentum cutoff is the 
same as replacing the one-body space L'^{M.^;C) by S^a- This procedure does not 
change the HFB energy 5hfb- But it changes the HFB evolution equation, because 



As usual Xa(^) is the characteristic function of the set A. 
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the constraint (|3.18p needs to be conserved for all times. The time-dependent 
equation in the cutoff model then becomes 



(3.19) 



«^rt = [PAFr,VA,Tt] 
at 

(7|t=o,a|t=o) = (70, "o) G A^HFB satisfying (I3.18p. 



where 

'Pa 



Va ■■= 



1-Pa 



The coupled system p.l3p can be rewritten in a similar fashion. 

Let us recall that the HFB dynamics introduced in the last section was obtained 
by restricting the (formal) unitary evolution of the Hamiltonian in Fock space to the 
submanifold of HFB states, by means of the Dirac-Frenkel principle. By introducing 
a cutoff in angular moment, we simply restrict ourselves to the time-evolution to 
the even smaller manifold of states belonging to fjA- We remark that, in numerical 
simulations, an angular momentum cutoff (as we introduced) will be in fact always 
imposed. 

The following main result of this paper now establishes blowup in HFB theory 
with an angular momentum cutoff. 

Theorem 4 (Blowup for HFB with angular momentum cutoff). Suppose w satis- 
fies assumption (12. 5p . Let (70,010) € ^hfb be a radially symmetric initial datum 
satisfying p.lSp and such that 

Tr|x|So + Tr(-A)7o <oo. 
Our conclusion is the following: // (70,0^0) has sufficiently negative energy, that is 
(3.20) SuFBho^cto) < -■^sup|w(r) +rw'{r)\^ ( (Tr7o)^ + Tr7o ) , 

then the unique maximal local-in-time solution (7f,af) of (j3.19p blows up infinite 
time, i. e., we have T < 00 and Tr(— A)"'^/^7t — s- cx) as t T~ . 

The proof of Theorem [3] is provided below in Section 14.31 

Remark 6 (Removing the angular momentum cutoff). It is not obvious how to get 
any information on the initial model from the study of the model with an angular 
momentum cutoff. However, it is easy to grasp from our proof in Section [4.31 that 
if Tr(|L|^+'^7t) stays bounded (or does not grow too fast), then blowup must occur 
under the same condition (|3.20p . 

Remark 7 (Case of the Laplacian). We mention, as a side remark, that in a 
nonrelativistic HFB model, where K is replaced by the Laplacian —A, and with 
the scaling critical two-body potential V{x~y) — — l/jx — yp, higher powers of the 
angular momentum are also a priori not conserved. However, in the nonrelativistic 
setting, one can use in the proof the observable M = jxp, instead of the operator 
AI — m? — Axi, as done below. Hence, for nonrelativistic systems, no 

angular momentum conservation is needed, and one easily sees that blowup always 
occurs under a negative energy condition, even for non-spherically symmetric initial 
data. 

Nevertheless, it should be mentioned that the Laplacian —A is qualitatively 
very different from the pseudo-relativistic operator A -I- m^, which models a 
relativistic system with a 'finite speed of propagation'. 
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4. Proofs 

4.1. Preliminaries. In the proofs of our main results, we shall need several aux- 
iliary lemmas, which we gather in this section for convenience. 

4.1.1. Spherical Harmonics Expansion. An important tool will be the expansion of 
the time-dependent states in terms of spherical harmonics. As we restrict ourselves 
to spherically symmetric states, our states will be uniformly distributed in each 
angular momentum sector and we will only encounter Legendre polynomials. 

We recall that the £-th Legendre polynomial Pe{t) is a real polynomial of degree 
£ defined on [—1; 1] which satisfies —l^Pe{t)^l and Pi{l) = 1, for all integers 
£ ^ 0. Its is linked to the usual spherical harmonics Y™ by the formula 

1 



(4.1) Pdco ■ co') = ^ Yr{cu)Yl-{oj'), 

7n——£ 

where eS'^ = {xeR^ : \x\ = 1}. 

We define for all £, £' the following function: 



F(;,t^'{r,r') = 27r ^ Pt{t)Pi,{t) V [yjr^ +r'^ - 2rr't^ dt 

(4.2) ^2tt [ Pi{t)Pij{t)( ^ ^» (\/r-2 + r'-i - ').rr't\ \ dt 

We recall that V{x) = — l/|a;| + ^dxl). With the usual abuse of notation, we will 
henceforth write T^(|a;|) as well as V{x), whichever is more convenient for the case 
at hand. 

The function Fi^ii will often appear when integrating out the angle variable in 
the interaction potential between the projections of our states in, respectively, the 
£-th. and the ^'-th angular momentum sector. Wc will need the following 

Lemma 1 (Estimates on Fiji). We have the following bounds: 

47r(l + \\r'w{r)\\r^,r, . ) 

(4.3) \F,j,{r,r')\^^ 

max(r, r') 

and 

(4.4) \r^ dr F,j.{r, r')| ^ C [l + £^ + {£'f) + ||(1 + ^')'+^^'(0|L=.[o.^) 
for some universal constant C and some constant depending only on e. 
Proof. Clearly we have 

IP t 'M<ro (^^\\ (Ml \ f' * 47r(l-f ||ru;(r)||^. ■ 

\Fij'(r,r)\i^2TT[l + \rw{r)\\^r^) - 



-1 ^r"^ -I- r'2 - 2rr't max(r, r') 

which is nothing else but Newton's Theorem. For the second bound (|4.4p . we treat 
w and the Newton interaction separately. Recall that 



+ r'2 - 2rr't ^„ max(r, r')'' 



-Pm(t). 



Differentiating explicitely we find 

/■I 

S$ C (l-f £2 _|_ (f 



r-3. r ^irn^dt 

J-1 Vr2 + r'2 - 2rr't 



^J2{l + m) I P,Pi,P, 
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where we have used that the integral J^-^ PiPe'Pm vanishes when m < \£ — £'\ and 



when m > i + £' , and that 



^ 2. For the term involving w, we write 



r'^dr j Pt{t)Pfj{t)w (Vr2 + r'2 - 2rr't^ dt 



dt 



s; \\{l + r'y+'w'{r) 



dt 



-1 (1 + r2 +r'2 - 2rr't) 



Then we explicitely compute 



dt 



(1 + r2 +r'2 - 2rr't)^+'' 2er' \{l + {r - r'fY (l + (r + r')2) 



1 



2rr' 



l+r2 + (r')' 



^ a 



2er' (l + (r-r')2)^(l + (r + r')2)'^ 

which ends the proof of Lemma [TJ 

The following is a simple consequence of Lemma [T] 
Corollary 1. Let p G L"'^(M'') be spherically symmetric function. Then we have 

(4.5) \V * p{x)\ < ^ ^^-^^ 

and 

(4.6) |V,|a:p(T/*p)(a;)| 

^ C (l + ||r«.(r)||^^(„_^) + a 11(1 + ''')'+^^'W|Lo.(o,oo) 
Proof. Integrating out the angle, we find 

V^p{r) = / {r'fdr' Fa^{r,r')p{r') 



□ 



hence the result is a simple application of Lemma [1] with £ — I' — Q and using that 
|V|x|2 (F * p)| 2\x\ \V * p{x)\ + |x|2 \VV * p{x)\. □ 

4.1.2. Pseudo-Relativistic Kinetic Energy. We now define the positive operator Ke^ 
acting on the space :— {[0 , oo) , r"^ dr) by 



(4.7) 



. 1 a 2 5 2 + 1) 



j,2 Qj, 



i.e. Ki is nothing but the restriction of \/ m2 — A to the i?-th angular momentum 
sector. We note that, for all k = —I, 



. . . , 1^, 



and hence 



v/m2 - ^Y^{lo,)u{\x\) = y/(w,)(X, w)(|a;|). 
We will need two important results involving the operator K^. The first one is the 
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Lemma 2 (Commutator Estimate). Suppose m > 0. We have, for all f € 



(4.8) [y^,/(x)J ,^,.(,3„^^I|V/II^- 

In particular, we have for all radial function f G W^'°°{M.'^) and all nonnegative 
integer £, 

(4.9) II [K,, fir)] ||^(^.([o,^),,.,,„ < C ||9./(r)||,.(„,^) . 

A proof of (|4.8p can be found in [5T]; see in particular the Corollary on page 309. 
The statement of this corollary does not give an effective bound on the norm of the 
commutator. However, the corollary is based on Theorem 3, on page 294 of [U], 
whose proof provides the effective control we need (a remark in this sense can be 
found in the paragraph 3.3.5, on page 305 of [21]). On the other hand, (|4.9p can be 
proven by applying (|4.8p in the subspace of L'^{R^) consisting of functions of the 
form {ujx)u{\x\) , which is stabilized by the operator [y/m'^ — A, /]. 

The second important result that will be very useful is the following 

Lemma 3 (Estimating Ki — Kfj). There exists a constant C such that one has, 
for all nonnegative integers 

(4.10) ll(^^-^£')i6(L^([0.oo),.^d.)) ^C{l + l + l')\£-t'\. 

The proof of Lemma [3] cannot be directly deduced from the literature, and we 
provide its proof in Section E] below. 

4.2. Hartree-Fock Theory: Proof of Theorem [H 

Step 1. Conservation of angular momentum. First of all, the proof of Theorem [2] 
makes use of the fact that p.2|) preserves both the spherical symmetry and the 
total angular momentum. We recall that L = —ix A V. 

Lemma 4 (Conservation of spherical symmetry and angular momentum). Let 70 € 
/Che be a spherically symmetric HF state. Then the unique maximal solution 74 to 
h2.1]) with 7|t=o — 7o is spherically symmetric for all times t G [0;T). 

If furthermore Tr(|Lp7o) < 00, then 7^ satisfies for all times [L,^t\ = and the 
total angular momentum is conserved: 

Tr |L|27t = Tr |L|27o for all t e [0,T). 

Note in particular Lemma[l]says that Tr|Lp7t < 00 for all t e [0, T). 

The proof of this lemma can be found in [T^] (it is trivial to see that the intro- 
duction of the regular and spherically symmetric potential w and the fact that 7 
may have infinite rank do not affect the proof of this lemma). The conservation of 
the squared angular momentum as expressed by Lemma |4] will be a key tool in our 
proof, as it was already in [T^] . 

For the rest of the proof, we introduce the non-negative self-adjoint operator 

3 



M :— Xi\/—A + Tn^ Xi 



1=1 



As in [inilll], the strategy of the proof consists in showing that the second derivative 
of the expectation TrM7( is negative, if the energy £HF(7t) = '?hf(7o) satisfies a 
certain condition. 
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Step 2. Estimate on direct and exchange terms. We will prove in this step that 
there exists a constant C (depending on w) such that 

(4.11) ^ Tr Af7t Tr A-ft + CTr(7o) Tr(l + \L\^)jo 
at 

with A = X ■ p + p ■ X. Here and henceforth, we use the notation p = — iV^ for 
the momentum operator and we use the notation C to denote any constant which 
only depends on w. The explicit dependence in w is very easy to derive, as it only 
occurs from the estimates of Lemma [TJ 

Remarks. One can check that Tr((l + |a;|'' + |pp)7t) < +cx) for t e [0, T), provided 
that initially Tr((l + ja;!"* + |pp)7o) < +oo holds. This fact may be used to show 
that all terms below are well-defined. Also, we note that 

Tr(M7t) ^ CTr(l + \x\'' + \p\'')^t) < oo 

and 

Tr(|^|7t) < CTr((l + \x\^ + \p\'')jt) < oo 
for aU times t £ [0,T). 

Remark 9. In order to have well-defined terms, we impose the higher regularity 
condition ||7o||'h2 < oo; see also Remark 1 iii) above. Note that a breakdown of the 
solution in || • 11-^= for some s ^ 1/2 implies breakdown of the solution in energy 
norm || • ||-^i/2; this claim follows easily by adapting the arguments in [T^. 

We start by computing 

(4.12) ^TrAf7t = TrA-ft - iTi [M, V * p^J jt + iTr [M, i?^J 7*, 
where we recall that 

R^{x,y) = V{\x - y\)-^{x,y) ^ '— + -f{x,y)w{\x - y\). 

\x-y\ 

To prove (|4.1ip . we have to bound the last two terms on the r.h.s. of the last equa- 
tion, arising, respectively, from the direct and the exchange term in the Hartree- 
Fock equation ()2.2p . We will estimate TrA7( later in Step 3. 

The appropriate estimate on the direct term is given in the following 

Lemma 5 (Estimating the direct term). We have for all spherically symmetric 
7 e /Chf 

(4.13) |Tr[Af,(T/*p^)]7| s;C(Tr(7))' 

where C only depends on w. 

For the exchange term, we have to face the problem that 74 can a priori be infinite 
rank, hence the method of [12] does not apply directly. However, an estimate similar 
to what was proved in |12j holds true, as expressed in the following 

Lemma 6 (Estimating the exchange term). We have for all spherically symmetric 
7 G /Chf such that Tr|Lp7 < 00 

(4.14) |Tr [Af, R^] 7K C Tr(7) Tr(l + \L\'')j 

where C only depends on w. 

Inserting the estimates (|4.13l) and (|4.14p in the expression (|4.12p . one gets the 
bound We now provide the 
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Proof of Lemma [31 Following the strategy of [TU] , we find 



(4.15) iTi [M,V * p^]j ^iTr ^/\p^ + ^ , \xf {V * p^) 

- 2ReTr 



V\p\'- 



which implies that 

|Tr [M,{V*p,)]jt\ ^ [\\[VW+^A^\'iV*p^)]\\+2\\\x\{V*p,)\\JjTT{j). 
Applying Lemma [2] to the first term in the parenthesis, we obtain 

|Tr [M,{V*p^)]j\ «C {\\W,\x\^V*p4oo + m{V*p^)\\oo)TT{j). 
Using the spherical symmetry of , (|4.13p is then a consequence of Corollary[TJ □ 

It rests to write the 
Proof of Lemma O First we write as before 



(4.16) iTi [M, R^] 7 = iTr ( Vbp + m^\x\^R^j 



iTr 



(i?^|x|VbP + ^'^'7) + 2ReTr (^R^x ■ -j=, 



We claim that 
(4.17) 

Tr (VbP + lkpi?^7) - Tr VbP + ^^7 



sc: CTr(7) Tr(l + \L\^h, 



and that 
(4.18) 



Tr RjX ■ 



V\p\'' 



^ C(Tr7)2. 



To prove (I4.17P and (|4.18p . we observe that, because of the spherical symmetry 
of 7, we can expand the kernel of 7 as 

(4.19) j{x, y)=J2 9ii\^l + ' 

where lOx = LUy = y/\y\ and Pg is the £-th Legendre polynomial whose formula 

was recalled in (|4.ip . 
We have 



Vbxp + lge{\x\, \y\)Pi{uJx ■ ujy) = Piioj^ ■ ujy) {Ke)rgii\x\, \y\) 

where Ki was defined in (j4.7p . The subscript x in and the subscript r = |a;| in 
{Ki)r indicate that these operators act on the x, respectively, r — \x\ variable. It 
follows that 
(4.20) 



Tr(A/bP+TO^|a;pi?77 



liy,x)y/\px\^ +m? {\x\^R^{x,y)) dydx 

= ^(2£+l)(2£' + l) / g,,{\yl\xm,{LOx-LOy)y. 
0, "'R" 



X ^J\px? ^-rn?^ {\xfgii\x\, |y|)Pf(wx • ^y)V(x - y)) dydx 
^(2£ + 1)(2£' + 1) (5£'(r,r'), {Ke^WFfj'ir, r')ge{r,r')) ^ 
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where Fi ii was defined before in dO]). Similarly to we obtain 

(4.21) Tr (^R^\x\^^\p\^+m^j) = ^(2^+ 1)(2/+ 1) {gej,F,j,r^ {Ki), g,) . 

Therefore we can rewrite the left hand side of (|4.17p in the form 

Tr [^/W^^\A^Ril) - Tr VIpP + l) 

= ^(2£ + l)(2f + 1) {ge,, {{K,,)^ r'F,,,, - r^F,,,, {K,\) 5^)^.^^. 
e,e' 

= ^(2^ + 1)(2/ +!)( (Ki - K^,)rr^F,^,,gi)^,^^, 

(4.22) + (.9,, [{K,)^yF,^,,] <?£')l2«l2 ) • 

Now using both Lemma [T] and Lemma [31 we get 
II [Kt - KtVFi^t II C(l + £ + 

This allows us to control the first term on the right side of (|4.22p . To control the 
second term, we remark that, from Lemma O we have 

\\[K,,\x\^F,,,,]\\ ^ \\dyFt,4^^^^^^.^, < C. 
Inserting these estimates in (I4.22p and noting that 

hdLl^Ll = hde^iLl) < ll5<?llsi(L?,) = TrL2g^ 

yields 

Tr(VbP + l|xpi?^7) - Tr VbP + "^'t) 

==; C^(2£ + 1)(2£' + 1) (1 + (1 + ^ + t)\t ~ e\) Tri2 gi, TrL^^gt 

1,1' 

^ C^(2£ + 1)(2/ + 1) (1 + ^2 _^ {I'f) TiLi gi TrL2g,, 

1,1' 

ii, C Tr7 Tr(l + |Lp)7 

and thus implies (|4.17|) . 

It remains to prove (|4.18p . Since x ■ p — rdr, with |a;| — r, we obtain 



Trii?. 



ge,Fe,e'rdr-^gi' 



< ^(2f +1)(2/ + 1) 

1,1' 

< ^(2f + 1)(2/ + 1) Wgde.^Ll) hAe.^Ll) \VFi,,'\\l^ \\dr{K,')-\^^,^^ 

1,1' 

< (Tr7)' 

as was claimed. Note that we have used that {dr)*dr = —r~^drr'^dr ^ K'^, for all 
£' ^ 0. Hence ||9r(^'^£')~'^ || ^ 1 as an operator acting on Lf. This ends the proof of 
Lemma [HI □ 
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Step 3. Estimating Ti (Ajt)- We recall that A = x ■ p + p ■ x. We will show in this 
step that 

(4.23) ^TrA7t < 2 £hf(7o) + (Tr7o)' sup \wir) + rw' {r)\_ . 

at 

In order to prove (|4.23p . we compute as before 

^ Tr Ajt = - I TrA[v/-A + + p^J - kR^^ , -ft] 

= - iTi [A, y/~A + vn? + n{V * p-yj - nR^Alt 
= 2 Tr ^ 7t 

With V{x) — — \x\ ^ + w{x), we obtain 

X ■ Wix) = - — r + X ■ Ww(x) — —V(x) + X ■ Vwfx) + w{x) . 
\x\ 



Since moreover p^ j \/p^ + "w? ^ -^/p^ + m^, we find 

4 Tr A-it 
dt ' 

<2fHF(70 

~ K. j dxdy {{x - y) ■ Vw{x - y) + w{x - y)) {-/tix, x)-/tiy,y) - |7f(x,y)p) 
< 2fHF(7t) + '«(Tr7t)^ sup|w(r) +rw'(r)|_ , 

where |/|_ ^ denotes the negative part of /. Here we used the fact that 
'yt{x,x)jt{y,y) ^ l7t(a^,?/)P since 7 0. Equation (|4.23p thus follows because 
fHF(7*) = ^hf(7o) and Tr(7t) = Tr(7o). 

Step 4- Conclusion of the proof of Theorem\^ It follows from Step 2 and Step 3, 
that 

(4.24) ^ TrM7t ^ f £:hf(7o) + ^ {^nof sup \w{r) + rw'{r)l 

+ t fTr(A7o) + C Tr(7o) Tr(l + \L\^)jo ] + TrMjo 



for all t < T (recall that T is the time of existence of the maximal local solution) 
and where C is a constant which only depends on w. From the assumptions on the 
initial datum 70, we find 

Tr(M7o) CTr(l + \x\^ + \p\^)-/o < 00 

and 

|Tr(A7o)| ^ \Tt{x ■ p + p ■ x)-fo\ ^ Tt{\x\^ + \p\^)-fo < . 

Hence it is clear that, if £hf(7o) + ('«/2)(Tr7o)^ sup^^Q \w{r) + rw'(r)|„ < 0, our 
estimate (I4.24p contradicts the positivity of the operator M, for t large enough. 
Therefore T < 00, and, from the blowup alternative (see Theorem [T]), it follows 
that ||7t||A'HF as < ^ T". □ 
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4.3. Hartree-Fock-Bogoliubov Theory: Proof of Theorem |4l The proof fol- 
lows the same lines as the one of the HF case and not all the details will be provided. 
We take as before M — Yl^i=i "vA^i^^^^s^i and note that M commutes with L, 
hence it also commutes with our cutoff projector Pa- The same calculation as in 
(Pl^ yields 



(4.25) ^Tr(M7t) = Tr([M, P^H^^PA\-ft) 



w{\x-y\) ]Pa®Pa 



at 



.\^-y\ 

However, since [M, Pa] — 0, and jt and at are easily seen to satisfy (I3.18P for all 
times, we may simply rewrite this as 



(4.26) ^Tr(M7t) = Tr A-/t - i Tr [M, V ^ p^,]-ft + iTr [M, i?^J jt 

1 

- ^n{at, [M^ + My,V{x - y)]at) ^2(^K6y 
With the help of Lemma [5] and we deduce 

\ 2 



and 



|Tr[M,y*p^j7i| s;C(Tr7t)^ 



Tr [M, i?^J 7t| CTr7t Tr(l + \L\'')-ft 



Now, we need an estimate on the pairing term, i. e. the last term of (|4.26p . This is 
the goal of the 

Lemma 7 (Estimating the pairing term). Let (7, a) G /Chfb be a radially symmet- 



HFB state such that Tr(|L|S+'^7) < 00 /< 



e > 0. The 



has 



(4.27) {a,[M, + My,Vix-y)]a) 



< C(Tr(l + 1^13)7) ^/^(Tr(l + 11^^!)''' (^^ TTi^) 



1/2 



where C only depends on w. 

Assuming that Lemma [7] holds true, we can conclude the proof of Theorem S) 
Combining estimates, we arrive at the bound 



^Tr(M7.) 



^ Ti-A-ft + CTr(70Tr(l + \L\^ht + CTr(l + |L|3)7,)^/'(Tr(l + \Lf+'ht) 

^TrA7i + c(l + A^) (Tr7t)' 

where we have used that Pa 7 = 7. 

The next step is to note that [Pa, ^] = 0, hence we have like in the HF case 



(4.28) -Tr^7t ^ 2 5HFB(7t,at) 

-K J dxdy {{x-y)-Vw{x-y)+w{x-y)) {■ytix,x)-/tiy,y) - \lt{x,y)\'^ + \atix,y)f) . 

We may use as before that '-ft{x,x)jt{y,y) — \7t{x,y)\'^ ^ and that ||at||^2(R6-) ^ 
Tr(7t) to infer 



dt 



Tr A-jt ^ 2 £hfb (7* , a* ) + k sup | w (r ) + rw' (r) | _ ( (Tr7t ) ^ + Tr7t ) 
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The end of the proof is then the same as in the HF case, using that Trjt = Tr7o 
for all times. 

It now rests to give the 
Proof of Lemma We write as usual 

(4.29) j{x, y)^J2 \y\)m + l)Pd^. ■ UJy) 

where ge{r,r') is the kernel of a self-adjoint operator ^ gi ^ 1 acting on radial 
functions in L^{W^;C), and 

(4.30) a(a;,y) = ^a,(|a;|,|y|)(2£ + l)P,(a;,-a;j,). 

This time the operator ai is antisymmetric, ai{r,r') — —ai{r',r), i.e. its kernel 
a^(|a;|, |y|) can be interpreted as a spherically symmetric, fermionic two-body wave- 
function. 

We also recall [5] that the condition ^ F ^ 1 gives aa* ^ 7(1 — 7). The 
same inequality must hold on each angular momentum sector, hence we obtain the 
following inequality for operators acting on L^: aia'^ ^ gg{l — gi). 

Now, arguing as for the exchange term, we calculate: 

(4.31) i (a, [M, + My, Vix ~ y)] a)^,^^,^ 

= 2i (a, + V(x - y)] a)^^^^^^ 

-mla,^=i^==-xV{x-y)a) 

\ VlP.P + l /i.(„6) 

For the second term on the right hand side of (I4.3ip we obtain 



(4.32) ( a, j§^rV{x - y)a 



V(2£ -f l)(2f + 1) ( ae, -^rF.^i.at 



where Fg^t' was defined before in (|4.2p . This yields to 

irai ^ ^(2£+ l)(2f + 1) ||adle,(i.) WWe^iLi) \VFiA\L-(o,oor ■ 
1,1' 

By our assumptions, we have aia*^ ^ gi{l — gi), hence 

Together with Lemma [1] we can further bound (|4.32p using the Cauchy-Schwarz 
inequality as 



\^M\^ |^^(2£+l)(Tri2^C^5,)i/2j <CTr(l + |L 
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Note the previous bound can be improved to s$ CTr (l + |L|2+'^) 7 for e > 0. 
Similarly as in the proof of Lemma [HI the first term in (14.31|) can be expressed as 



(4.33) ^(2^ + l){2l' + 1) {a,,, ( (if,,). " r^F,,^, {K,)^)a,,) 

= ^(2^+ l)(2f + 1) {a,,, {K,, - Ke)yFe>jai)^,^^, 

+ ^(2£ + 1)(2£' + 1) (a,,, [(Ki)^ , r^F,,,e]ai)^, 

1,1' 

Using Lemmas [l] [5] and [3] like for the exchange term, we obtain 
^(2^+ 1)(2/ + 1) (a,,, {Ki,\x\''Fti - |a;pF,,,,if,)a,,)^,^ 

^ C J2i^£ + 1){2£' + 1)(1 + {£ + e')\£ - £'\) ii-r „e 

e,e' 

^c(^{2£+l){l + e^)TTL2{gey^'^ |^X^(2/ + l)Tri.(g,0'/'^ 

^ C(Tr(l + |Ln7)^/^(Tr(l + iL^^h)'^' (^^ TtW) 
This ends the proof of Lemma [T] □ 
Appendix A. Proof of Lemma [3] 



1/2 



We have, using the well-known [15^ integral formula for -^/—A + — m and 
assuming that w, v are smooth enough. 



= ((v/m2 - A - m)Y,\cj,)ui\x\), r/(c^,)i>(|a:|) 
1 ff IC2i\x-y\) 



^^'JU k-yP -{Yt{Lu,H\x\)~Y,H^y)ui\y\))x 
X (r/(^.)«(|x|) - r/(^,)z;(|y|))dxdy 
where IC2 is a Bessel function. Averaging over k and using that 

^ ^^r/(c.')-PK^-c^') 

fc=-£ 

where is the -^-th Legendre polynomial, we get 

(A.l) {{Ki~Ki,) ru,w)i2([o,oo),r2dr) 

" "2^ y ije [x-yP (Pd^x ■ Wy) - Pe'{uj^-ujy))dxdy. 

Note that we have used that for all £, Pe{ujx -^^x) — Pe{^) = 1- Using that |a;p/C2(a;) 
is bounded, this yields 

(A.2) {{Ki -Ki,)r u, v) ^^(^Q^^y^^ar) 

1 [[ \x\\u(\x\)\\v{\y\)\\Pi{uJx ■ i^y) ~ Pl'{i^x ■ t^y)\ 
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Passing first to spherical and then to polar coordinates gives 

(r^ _l_ g2 _ 2rsa) 



-1 JQ 



4 1 da I t^dt I ^'^.cos3.sin^.l-(^^°^^)l "^(^^-^^l 



-1 Jq Jq (1 — 2 sin 61 cos 

Using the Cauchy-Schwarz inequality for the t integration we get 

. f' , ,.cos3/2 0sini/2 0|Pf(a)-Pf,(a)| 

J^i Jo [1 ~ 2 sm cos ua) 

Let us recall that ^^(l) = 1 for all £ and introduce 

\P,{a)-Pt{a)\ 



Ci^ii := sup 

such that 



qg[-i,i] 1-a 



^ II ,\ ^cos3/2 0sini/2 0(l-a) 



^ ■ '■ " J _i Jq (1 — 2sin6'cos6'a) 

which is finite as seen by computing 

1 1_« (l + u)log(l±H)-2u 

da 5- = 2 - log 1 - u). 

Hence it remains to show that 

(A.3) C^i^ i^C{l + e + e')\£~£'\. 

The Legendre polynomials satisfy the relation 

(1 - x^)P;,{x) = -nxPn{x) + nPn-i{x) 
which may be written 

1 — X n 
Assuming for convenience I' > £ and summing over n = ^ + 1, . . . , i?', we obtain 

P,,{x)-P,{x) _ 4^ (^^^^ , M^) 



1-x 

n=e+i 



E (Pn{x)-{1 



We have the formula 

n{n + l) 
l-a;2 

which shows that {1 + x)\P|^{x)\ ^ n{n + 1) . Hence we obtain the bound 



n=t+l 

which completes the proof of Lemma [2] □ 
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